Reversed propagation dynamics of Laguerre- Gaussian beams in left-handed materials 
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On the basis of angular spectrum representation, the reversed propagation dynamics of Laguerre- 
Gaussian beam in left-handed materials (LHMs) is presented. We show that negative phase velocity 
gives rise to a reversed screw of wave- front, and ultimately leads to a reversed rotation of optical 
vortex. Furthermore, negative Gouy-phase shift causes an inverse spiral of Poynting vector. It is 
found that the Laguerre-Gaussian beam in LIfMs will present the same propagation characteristics 
as the counterpart with opposite topological charges in regular right-handed materials (RHMs). The 
momentum conservation theorem insures that the tangential component of the wave momentum at 
the RHM-LHM boundary is conserved. It is shown that although the linear momentum reverses its 
direction, the angular momentum remains unchanged. 
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I. INTRODUCTION 



Almost 40 years ago, Russian scientist Victor Vese- 
lago proposed that a material with electric permittivity 
e < and magnetic permeability /i < 0, would reverse 
all known optical properties He termed these media 
as left-handed materials (LHMs) since the wave vector 
k, forms a left-handed triplet with the vectors E and H. 
That is, phase velocity and the Poynting vector are an- 
tiparallel, which consequently results in counter-intuitive 
phenomena such as reversals of the conventional Doppler 
shift and Chcrcnkov radiation as well as reversed refrac- 
tion. Veselago pointed out that electromagnetic waves 
incident on a planar interface between a regular right- 
handed material(RHM) and a LHM will undergo negative 
refraction. Hence a LHM planar slab can act as a lens 
and focus waves from a point source. Recently, Pendry 
extended Veslago's analysis and further predicted that a 
LHM slab can amplify evanescent waves and thus behaves 
like a perfect lens Pendry proposed that the ampli- 
tudes of evanescent waves from a near-field object could 
be restored at its image. Therefore, the spatial resolu- 
tion of the superlens can overcome the diffraction limit 
of conventional imaging systems and reach the sub- wave 
length scale. The physical realization of such a LHM was 
demonstrated only recentlyfor a novel class of engineered 
composite materials 0, 13,11) 0- 

After the first experimental observation of negative re- 
fraction, intriguing and counterintuitive phenomenon in 
LHMs, such as amplification of evanescent waves 0,01 j 
unusual photon tunneling [1, [^, and negative Goos- 
Hanchcn shift (Tol . IT]| have attracted much attention. 
Here we want to explore the reversed propagation dynam- 
ics of Laguerre-Gaussian beams in LHMs. The propaga- 
tion of Laguerre-Gaussian beam has been investigated in 
conventional RHMs [12, llj, . Such beams have a 



phase dislocation on the beam axis that in related liter- 
ature is sometimes referred to as an optical vortex [l^ . 
For a general Laguerre-Gaussian beam the Poynting vec- 
tor has an azimuthal component. This means that there 
is an energy fiow along the circumference of the beam as 
it propagates, giving rise to an orbital angular momen- 
tum [l^ . It is found that the spiral of the Poynting vec- 
tor of a Laguerre-Gaussian beam is proportional to the 
Gouy-phase shift [13, [Hi. It is known that an electro- 
magnetic beam propagating through a focus experiences 
an additional n phase shift with respect to a plane wave. 
This phase anomaly was discovered by Gouy in 1890 and 
has since been referred to as the Gouy-phase shift p^ . 
Because of the negative index, however, we can expect 
an reversed Gouy-phase shift in LHMs ^2u\ . Hence it will 
be interesting for us to describe in detail how the Poynt- 
ing vector evolves as it propagates and how the reversed 
Gouy-phase shift affects its spiral in LHMs. 
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In this work, we will reveal reversed propagation dy- 
namics of Laguerre-Gaussian beam in LHMs, such as 
inverse screw of wave- front, inverse spiral of Poynting 
vector, and inverse rotation of vortex field. First, start- 
ing from the representation of plane-wave angular spec- 
trum, we obtain the analytical description for a Laguerre- 
Gaussian beam propagating in LHMs. Our formalism 
permits us to introduce the reversed Gouy-phase shift to 
describe the wave propagation. Next, we will recover how 
the wave- front and Poynting vector evolves, and how the 
reversed Gouy-phase shift affects their propagation be- 
havior. Then, we attempt to investigate how the nega- 
tive index influences the linear momentum and angular 
momentum of Laguerre-Gaussian beams. Finally, we will 
explore how the negative index gives rise to the reversed 
rotation of the vortex field. For a comparison, the corre- 
sponding propagation characteristics in RHMs will also 
been discussed. 
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II. THE PARAXIAL PROPAGATION OF A 
LAGUERRE-GAUSSIAN BEAM 

To investigate the propagation dynamics of a Laguerre- 
Gaussian beam in LHMs, we use the MaxweU's equations 
to determine the field distribution both inside and out- 
side the LHM. We consider a monochromatic electromag- 
netic field E(r,t) = i?e[E(r) exp(— zwi)] and B(r, i) = 
i?e[B(r) exp(— icji)] of angular frequency oj propagating 
from the RHM to the LHM. The field can be described 
by Maxwell's equations 



V X E 

V X H 



(9B 

~~dt 
dD 



, B = fiQti ■ H, 

D-£o£-E. (1) 



From the Maxwell's equations, we can easily find that 
the wave propagation is only permitted in the medium 
with e, /i > or e, /i < 0. In the former case, E, H and k 
form a right-handed triplet, while in the latter case, E, 
H and k form a left-handed triplet. 

We introduce the Lorentz-gauge vector potential to 
describe the propagation characteristics of Laguerre- 
Gaussian in RHMs and LHMs. The vector potential of 
the beam propagating in the +z direction can be written 
in the form 

A = Ao{aej: + f3ey)up,i{r) exp{ikz — icut), (2) 

where Aq is a complex amplitude, e^^ and are unit 
vectors, k = tij^ luj/c, c is the speed of light in vacuum, 
n-R = \/£rtP-R. and = —^/eiJil are the refractive index 
of RHM and LHM, respectively The coefficients a 
and /3 satisfying a — i{af3* — are the polarization 

operators with a — ±1 for left-handed and right-handed 
circularly polarized light. 

When the field distribution is specified at a boundary 
surface or a transverse plane, one can obtain a unique 
solution of the electric field of the wave propagating in 
the +z direction. Here, we assume that the transverse 
electric field at the z = plane is given by a Laguerre- 
Gaussian function as follows: 



u{r,ip,0) 



Wq 



^/2r 



p 



2r' 



exp 



— ilip 



•(3) 



A Laguerre-Gaussian beam has two mode indices to fully 
describe the mode: / and p. A given mode will have I 
complete cycles of phase 27r upon going around the mode 
circumference, so that I is known as the azimuthal in- 
dex. The index p gives the number p + 1 of radial nodes. 
Laguerre-Gaussian light beams are well known to possess 
orbital angular momentum due to an exp[zZ(/3] phase term, 
where is the azimuthal phase. This obital angular mo- 
mentum Ih is distinct from the spin angular momentum 
due to the polarization state of the light [IG]. 

From the point of view of Fourier optics, we know that 
if the Fourier component at the z ~ plane represents 
the angular spectrum that the transverse component of 



the wave propagating in the half space z > should have. 
Then, the field in the region z > can be expressed by 
an integral of the plane wave components associated with 
the angular spectrum given at the z > plane '^I'l. The 
angular spectrum is related to the boundary distribution 
of the field by means of the relation 



u(fc) 



drrJi {kr)u(r, (p, 0), 



(4) 



where J/ is the first kind of Bcsscl function with order I. 
The two-dimensional Fourier transformations of Eq.Q 
can be easily obtained from an integration table ^^2] . In 
fact, after the field on the plane z = is known, Eq. ([3]) 
together with Eq. ^ provides the expression of the field 
in the space z > 0, which yields 



u(r,ip,z) 



dkk exp 



Ji{kr)u{k). (5) 



which is a standard two-dimensional Fourier trans- 
form The field ii(r, z) is the slowly varying envelope 
amplitude which satisfies the paraxial wave equation 



1 



d 



dz 2n[irko 



vi 



u{r,z) = 0, 



(6) 



where Vj_ = dxGx + dyBy. From Eq. ^ we can find that 
the field of paraxial beam in LHMs can be written in 
the similar way to that in RHMs, while the sign of the 
refractive index is negative. 

The gauge condition on the vector and scalar potentials 
takes the form = (i/ k)W ■ A. The electric and magnetic 
fields are obtained from the potentials as 



E(r,t) 



dA 



du du^ 
a— + 13— le^ 

oy. 



dx 



iuj{aex + Pey)u 

exp(ifcz — iiut), (7) 



B(r,t) 



V X A = Aq 



du du 

f^a 

ox oy 



— ik{j3ex — aey)u 

exp{ikz — iujt). (8) 



These field expressions neglect terms in each component 
that are smaller than those retained in accordance with 
the paraxial approximation [23!|. The z components are 
smaller than the x and y components by a factor of or- 
der l/kwQ. It is readily verified that the fields satisfy 
Maxwell's equations. Note that the Cartesian derivatives 
can be converted to polar r and derivatives in the usual 
way. 

To be uniform throughout the following analysis, we in- 
troduce different coordinate transformations z*(i = 1,2) 
in the RHM and the LHM, respectively. First we want 
to explore the field in the RHM. Without any loss of 
generality, we assume that the input waist locates at the 
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object plane z = —a and = z + a. The field in the 
RHM can be written as 



a 



pi 



w{zl) 



X exp 



^/2r 

^2(2*) 

injikozl 



p 



y/2r 

inRk^r'^zl 



exp 



exp[— 



X exp[— i(2p 



R{zl) 

1) arctan(0j/zi^)], (9) 



wiz 



{zI)^wo^1 + {zI/zr)\ R{zl)^zl + 



(10) 



Here Cpi is the normalization constant, L^^[2r'^ / w\{z*y\ 
is a generalized Laguerre polynomial, zr = nRkQwl/2 is 
the Rayleigh length, is the beam size and R{zl) 

the radius of curvature of the wave front. The last term 
in Eq. ^ denotes the Gouy phase which is given by 
$1 = -(2p+ 1^1 + l)arctan(zi*/zfl.). 

We are now in a position to calculate the field in LHM. 
In fact, the field in the RHM-LHM boundary can be 
easily obtained from Eq. ^ by choosing z — 0. The 
plane-wave spectrum of the Laguerre-Gaussian beam can 
be obtained by performing the two-dimensional Fourier 
transform in Eq. Q. After the plane- wave spectrum on 
the plane z = is known, Eq. ([5]) provides the expression 
of the field in the space z > 0. For simplicity, we assume 
that the wave propagates through the boundary without 
refiection, the field in the LHM can be written as 



Co 



w{z*) 



X exp 



V2r 

in^koz^ 



4' 



V2r 



exp 



exp[ 



X exp[— i(2p 



R{z*) 

1) arctan(2:2/-2L)] 



wHz*) 
ilip] 

(11) 



w{z;) = wo Ji + i4/^, Ri4)^4 + ^- (12) 

^2 

Here Z2 = z — (1 — nL/riRja and zl = nLkoWQ/2 is the 
Rayleigh length in LHM. The beam size ^(^1) ^^'^ the 
radius of curvature R{z2) are given by Eq. (fT2)l . The 
Gouy-phase shift in LHM is given by = — {2p -I- |^| -I- 
1) arctan(z2/zL). Because of the negative index, the re- 
versed Gouy-phase shift should be introduced. A more 
intuitive interpretation of the reversed Gouy phase can 
be given in terms of a geometrical quantum effect (23 | 
or the uncertainty principle [1^. As can be seen in the 
following section, the inverse Gouy-phase shift will give 
rise to an inverse spiral of Poynting vector. 

For a Laguerre-Gaussian beam with I ^ 0, the on-axis 
phase form exp[il(p] results in that the surfaces of wave- 
front have helical form. Specifically, I refers to the num- 
ber of complete cycles of phase 27r upon going around 
the beam circumference. Now let us to study the screw 
of the wave front. Here the sense of the positive angles 




RHM 



FIG. 1: (Color online) The helical wave front for Laguerre- 
Gaussian beam with I = 1 result from an azimuthal phase 
structure of exp[— iiys]. In the LHM, the phase velocity Vp 
reverses its direction. The wave-fronts exhibit anti-clockwise 
screw in the RHM, while present clockwise screw in the LHM. 



is chosen as anticlockwise, while negative angles are con- 
sidered in the clockwise direction. In the regular RHM, 
the constant wavefront satisfies 



nRkoz'l 



-iuRkpr^ 
R{zl) 



lip -\- <^Q = const. 



(13) 



The schematic view of the wave front is a three- 
dimensional screw surface of (rcos(p, rsin(/j, z). The 
plotting range of r is from to 5wo with the interval of 
Ar = O.Swq and that of URk^z is from — 47r to with 
the interval of nj^fcoAz = O.Itt. The wavefront structure 
exhibits a anticlockwise-screw type with a pitch of Xq/ur 
along the +z axis. 

Next we explore the screwing fashion of wave front in 
the LHM. The constant wavefront satisfies 



riLkozX + 



lip + = const. 



(14) 



The plotting range of r is from to Swq with the in- 
terval of Ar = Q.^wq and that of ULk^z is from to 
—Ait with the interval of n^^fcoAz = —O.Itt. The wave- 
front structure is a clockwise-screw type with a pitch of 
•^o/I^-lI along the -\-z axis. Figure [2] shows a typical form 
of a helical wavefront structure from the RHM to the 
LHM. At the RHM-LHM interface, the wave front will re- 
verse its screwing fashion. It is intriguing to observe that 
the wave-front of Laguerre-Gaussian beam with I = \ m 
LHMs will exhibit the same skewing fashion as the coun- 
terpart with I = —1 h\ RHMs. As can be seen in the 
following section, the inverse screw of wavefront will re- 
sult in an inverse rotation of optical vortex. 
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III. POYNTING VECTOR AND ANGULAR 
MOMENTUM 

The propagation characteristics of electromagnetic 
fields are closely linked to their local energy flow, which 
is usually discussed by use of the Poynting vector. There 
has been considerable interest in orbital angular momen- 
tum of Laguerre-Gaussian beams relating to Poynt- 
ing vector in free space. The Poynting vector has a mag- 
nitude of energy per second per unit area and a direction 
which represents the energy flow at any point in the field. 
The time average Poynting vector, S can be written as 



1. 



-RefE X H* 



(15) 



The spiral of the Poynting vector in free space or regular 
RHMs has been discussed extensively [H, [13, H 0, [13, 

m. 

Now a question arise: what happens in LHMs with 
simultaneously negative permeability and permittivity? 
The potential interests encourage us to derive a general 
expression to describe the Poynting vector in RHMs and 
LHMs. Substituting the expression of Eqs. ([7]) and ([8]) 
into Eq. (fT5|) we find 



Sip 



1 Lokr , 



^flQ R 

1 



1 



(16) 



Here the component Sr, relates to the spread of the beam 
as it propagates. The azimuthal component S^p describes 
the energy flow that circulates around the propagating 
axis. The presence of this flow is due to the existence 
of the longitudinal components and of the field. 
The first term of the azimuthal component depends on /, 
where Ih has been identified as the orbital angular mo- 
mentum per photon |l^. Its second term relates to the 
contribution of polarization and intensity gradient. The 
contribution of circular polarization will lead to a spin 
angular momentum of the beam. The axial component 
Sz describes the energy flow that propagates along the 
+z axis. 

Next, we attempt to explore the angular momentum 
in LHMs. Since the dispersion cannot be ignored in a 
causal system with a negative index of refraction. The 
momentum conservation theorem should be derived from 
the Maxwell equations and the Lorentz force and is given 

by di 



V • T- 



dG 
'dt 



(17) 



where G is the momentum density vector and F is the 
force density. The momentum flow T also referred as the 
Maxwell stress tensor. It is well-known that a material 
contribution to the energy density accompanies the prop- 
agation of electromagnetic energy in dispersive materials. 



Analogously, there exists a corresponding material con- 
tribution to the wave momentum. Thus the momentum 
conservation equation for the electromagnetic wave can 
be written in the form 1301 



G 
T 
F 



.,.E.H.4(||EP + ||HF 



^Re 
2 

iRe[(D • E* -(- B • H*)/ - (DE* + BH*)], 
iRe[peE* + J X B* -f p„H* + M x D*] . (18) 



Here the momentum density G contains the Minkowski 
momentum Gm = D x B plus material dispersion terms. 
The tensor J is 3 x 3 identity matrix. The electric and 
magnetic polarization vectors are give by Pg = eo(e— 1)E 
and P„i = — /io(/i — 1)H, respectively. Bound electric 
current J — dV^/dt and bound electric charge pe — V-Pg 
have been accounted. Similarly, bound magnetic current 
M = dPm/dt and bound magnetic charge Pm. = ^ ■ Pm 
should be introduced to describe the angular momentum 
flow in LHMs. 

Now we want to enquire: how to determine directions 
of the momentum density and the momentum flow? It 
is well known that the time-domain energy density in 
a frequency nondispersive medium are defined as = 
i [D • E + E • H] . Obviously, the energy density in LHMs 
would be negative if the permittivity and permeability 
were negative. Hence, the energy density in a frequency 
dispersive medium is defined as [31, 32] 



d{eLj) 
du! 



El 



d{puj) 
du! 



H 



(19) 



In principle, the energy density can be decomposed 
into electric and magnetic parts. The positive elec- 
tric and magnetic energy requires d{eLj)/duj > and 
d{pLu)/duj > 0. Subsequent calculations of Eq. (fT5|) show 
that both the momentum density G and the momen- 
tum flow T in LHMs are antiparallel to the power flow 
S — iRe[E X H*]. Hence we conclude that the linear an- 
gular momentum flux will reverse its direction in LHMs. 
Note that Eqs. (fTS]) and (fT9)) are valid only for lossless 
media, and its application to lossy media produces un- 
physical phenomena such as a negative energy in LHMs. 
In a lossy and dispersive LHM, the momentum flow of a 
monochromatic wave is opposite to the power flow direc- 
tion. However, the momentum density may be parallel 
or antiparallel to the power flow [33] . 

The cross product of this momentum density with the 
radius vector r yields an angular momentum flow. The 
angular momentum flow in the z direction depends upon 
the component of G^, such that 



J. = r X 



(20) 



Conservation of momentum at a material boundary en- 
sures that the tangential component of the wave momen- 
tum is conserved [30 . 33J . Hence the angular momentum 
flow still remains unchanged in the LHM. We can predict 
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FIG. 2: (Color online) Ray optical model of Laguerre- 
Gaussian beam in the RHM and the LHM. The rays (green 
arrows) lies on a single-sheeted hyperboloid surface. Note 
that the arrows indicate the direction of the Poynting vec- 
tors. When Rays travel from a RHM to a LHM, the negative 
refraction results in the reversed screw. The rays exhibit anti- 
clockwise screw in the RHM, while he rays present clockwise 
screw in the LHM. 



theoretically that the orbital angular momentum for per 
phonon still remain Ih. In order to accurately describe 
the angular momentum flow, it is necessary to include 
material dispersion and losses. Thus a certain disper- 
sion relation, such as Lorentz medium model, should be 
involved. 

It has been shown that the trajectory at peak inten- 
sity becomes a straight line skewed with respect to the 
beam axis in RHMs [34| . To obtain a better physical pic- 
ture of the straight trajectory in LHMs, the ray optical 
models of Laguerre-Gaussian beam are plotted in Fig. [2l 
Within a ray optical picture, the angular spectrum may 
be represented by skew rays in the optical beam. The 
screwing behavior of rays can be deduced from Eq. (jl6[) 
in which we see that each ray having an azimuthal an- 
gle = l/{nB.,Lkor) and a polar angel rj = r/R with 
respect to the beam axis. Hence all rays lies on a single- 
sheeted hyperboloid surface. Note that wave-vector and 
the Poynting vector are parallel in the RHM and antipar- 
allel in the LHM. Energy conservation requires that the 
z component of Poynting vector must propagates away 
from the interface. Both wave-vector and energy flow in- 
cident on a planar interface between a RHM and a LHM 
will undergo negative refraction. Therefore all rays will 
reverse its screwing fashion in the LHM (see Fig. [2|) . 

To explore the reversed propagation dynamics in the 
LHM, we need to move beyond ray tracing. The ray 
optical models neglect diffraction and thus could not be 
used to predict precisely the spiral of Poynting vector. 
To include diffraction, we had to use a more accurate de- 
scription of the electromagnetic field. It is well known 
that the trajectory of the Poynting vector is described 



by the spiral curve The relative value of the com- 

ponents determines the trajectory of the Poynting vec- 
tor. The spiral angle of the Poynting vector is given by 
9 = 6o + Kz, where k = dO/dz = S^p/{rSz) is the rate of 
rotation. The period of the trajectory along the z axis 
is lix I K. For a general Laguerre-Gaussian beam, the rate 
of the azimuthal rotation is relate to the distance given 

by 

dz n_R,Lfcor2 riR^Lkor^ nR,Lkow'^{z) 

^ 4a L^;;ll^[2rywHz)] ^^^^ 
riR.LkQW^iz) L|!'[2rVu;2(z)] 

Note that the first term is polarization independent. 
While the last three terms depends on the polarization. 
For modes with p = 0, the final term is always zero. 

For a single-ringed Laguerre-Gaussian beam p — Q but 
Z 7^ 0, the scaled radius of the peak intensity is given by 
Tmax = vW/2w(^)- We find that, for all values of I and 
a, the rotation angle is given by 

I z 
Ornax = TTT arctau . (22) 

I' I ^B.,L 

Figure [3] shows the vector fields illustrating the spiral an- 
gle of the Poynting vector. The Poynting vector exhibits 
anticlockwise spiral in the RHM as depicted in Fig.[3l[a), 
while presents clockwise spiral in the LHM as plotted in 
Fig. ^h) . The theoretical analysis and numerical calcu- 
lations presented here coincide with experimental obser- 
vations [sH]. For / ^ 0, Laguerre-Gaussian beam have 
annular intensity profiles and as rmax is typically much 
greater than optical wavelength A, the skew angle is ex- 
pected to be very small. 

It can clearly be seen that the sense of spiral (clock- 
wise or anticlockwise) depends on the signs of I and the 
Rayleigh length. However, the amount of rotation of the 
poynting vector is independent of the magnitude of I. It 
is interesting to note that the Laguerre-Gaussian beam in 
LHMs will present the same fashion of spiral as the coun- 
terpart with opposite topological charges in RHMs. This 
is consist with the ray optical model describing the spiral 
of the Poynting vector. Equation implies that the 
absolute rotation for a single-ringed Laguerre-Gaussian 
beam at z — zrx is 7r/4, regardless of I. When the far- 
field pattern of Poynting vector is calculated it is found 
that the Poynting vector is spiraled by 7r/2. 

When I — and p = 0, the Laguerre-Gaussian is iden- 
tical to the fundamental Gaussian beam, the spiral of the 
Poynting vector arises from the effect of circular polar- 
ization [18[. The azimuthal rotations in the RHM and 
the LHM are relate to the distance can be obtained as 
9 — (7 arctan(z/zfl) and 9 = cr arctan(z/zi), respectively. 
It is intriguing to note that the right-handed circularly 
polarized beam in LHMs will present the same fashion 
of spiral as the left-handed circularly polarized beam in 
RHMs. Hence we can describe quantitatively the amount 
of spiral of the Poynting vector, from which we can de- 
termine whether a material is a LHM or a RHM. 
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(a) 



FIG. 3: (Color online) Numerically computed field intensity 
distribution and transversal components of Poynting vector 
(green arrows) for / = 3 Laguerre- Gaussian beam. (a)Tlie 
Poynting vector exhibits anticlockwise spiral in RHM Zi = 
nnkoWQ /2. (b) The Poynting vector presents clockwise spiral 
in LHM Z2 — \nL\kQw'^/2. For the purpose of comparison, we 
have chosen til = —riR. 
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FIG. 4: (Color online) Interfering a Gaussian beam and a 
Laguerre- Gaussian beam of azimuthal index I — 3 will pro- 
duce the vortex field with three spiral arms, (a) The field 
distribution at plane of Zi = nijfcoWo/2. (b) The field distri- 
bution at plane of Z2 = |nL|fcoWo/2. When the vortex field 
enters the LHM, the rotation changes its fashion. 



Interfering the Laguerre-Gaussian beam with a fun- 
damental Gaussian beam will transform the azimuthal 
phase variation of the pattern into an azimuthal inten- 
sity variation. Hence the helical phase ultimately results 
in an vortex field with |Z| spiral arms [s^, [s^, H^. The 
intriguing properties strongly motivate us to explore the 
vortex field propagation in LHMs. A simulation of the 
vortex field rotation produced by interfering a fundamen- 
tal Gaussian beam and an Laguerre-Gaussian beam of 
azimuthal index Z = 3 is shown in Fig. [H It can be 
seen that the vortex has three spiral arms, which is a 
result of the mismatch between the wave-fronts of the 
Laguerre-Gaussian beam and the Gaussian beam. The 
reversed screwing wave-fronts will directly cause an in- 
verse rotation of the vortex filed in the LHM. The vortex 
will always has a spiral shape unless the wavefronts of 
the two beams have the same curvature. For example, 
the vortex field at a focusing waist will exhibit |Z| intense 



spots. After the vortex propagating through the focusing 
waist, the spiral arms will change their shapes. 

Now let us consider how to modulate the rotation 
of the vortex field. As we change the path length of 
the Gaussian beam, the spiral arms will rotate around 
the propagating axis. This is analogous to altering the 
phase different between the Laguerre-Gaussian beam and 
the Gaussian beam \3§\. The spiral arms repeat every 
A/ItilI in the LHM, but only rotate fully after propa- 
gating |Z|A/|nL|. A path length change in the Gaussian 
beam of 3A/|ni| will cause the pattern to rotate through 
27r and — 27r in RHM and LHM, respectively. The vortex 
presents an anticlockwise rotation in the RHM, while ex- 
hibits a clockwise rotation in the LHM. Once the vortex 
field enter the LHM, it will reverse its rotation fashion. 
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IV. CONCLUSIONS 

In conclusion, we have investigated the reversed prop- 
agation dynamics of Lagucrrc-Ganssian beam in LHMs. 
We have introduced the concepts of negative Gouy-phase 
shift to describe the propagation of Laguerre-Gaussian 
beam in LHMs. The negative phase velocity and nega- 
tive Gouy-phase shift caused inverse screw of wave- fronts, 
reversed spiral of Poynting vector, and inverse rotation 
of vortex field. At a RHM-LHM interface, direct calcula- 
tion of Maxwell's equations dictates the wave- vector and 
energy flow undergoes negative refraction. Consequently, 
inside the LHMs, the screw of wave- fronts, the spiral of 
Poynting vector, and the rotation of vortex will reverse 
their direction. We have shown that the Poynting vec- 
tor of Laguerre-Gaussian beam in LHMs will present the 
same fashion of spiral as the counterpart with opposite 
topological charges in RHMs. Conservation of momen- 
tum at the boundary ensure that the tangential compo- 



nent of the wave momentum is conserved. We have found 
that although the linear momentum reverses its direction, 
the angular momentum still remains unchanged. Since 
the photons in Laguerre-Gaussian beam possess angular 
momentum, the reversed propagation dynamics may of- 
fer new fundamental insights into the nature of LHMs. 
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